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Abstract 

We present a self-consistent analytic theory of the intra-layer and inter-layer pair correlation 
functions in electron-electron and electron-hole fluid bilayer systems. Our approach involves 
the solution of a zero-energy scattering Schrodinger equation with an effective potential which 
includes a Bose-like term from Jastrow-Feenberg correlations and a Fermi term from kinetic 
energy and exchange, tailored to yield the Hartree-Fock limit at high carrier density. The theory 
is also shown to satisfy the plasmon sum rule and the charge neutrality condition. We obtain 
good agreement with the available Quantum Diffusion Monte Carlo (DMC) data in symmetric 
bilayers over a wide range of carrier density and layer spacing, and stress the role of three- 
body correlation terms in yielding fully quantitative agreement at strong coupling. Signals of 
impending transitions to density-modulated states at strong coupling and low layer spacing 
appear in the calculated pair correlations through inter-layer in-phase oscillations for electron- 
hole bilayers and out-of-phase oscillations for electron-electron bilayers, in agreement with the 
DMC findings. 
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1. INTRODUCTION 

Advances in nanoscale growth technology have for a number of years allowed the 
fabrication of materials formed by alternating layers of two or more semiconductors. 
Electron or hole layers trapped at the various interfaces have made available Fermi- 
liquid systems having a quasi-two-dimensional (2D) low-energy dynamics under good 
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control of carrier density and inter-layer spacing [IJ. It was soon realized that systems 
of spacially separated carriers present some advantages over conventional samples in 
which the carriers occupy the same region. Thus in electron-electron bilayer systems 
(EEB) evidence has been reported for novel fractional quantum Hall states [2] and for 
a bilayer quantum Wigner solid [3|4j . In electron- hole bilayer systems (EHB), on the 
other hand, the Coulomb attractions between the two kinds of carriers can bring about 
inter-layer coherence and hence again crystalline order as well as exciton condensation 
and superfluidity [5]. 

Many of the Coulomb interaction effects in these systems can be understood with the 
help of the homogeneous 2D electron gas model [6] . A central role in the theory is played 
by the intra-layer and inter-layer pair correlations as expressed by the instantaneous pair 
distribution functions gafsif), where a Greek suffix a = 1 or 2 stands for carriers in 
one of the two layers. These functions are defined so that the quantity 2i:rnf3ga[3{r)dr 
gives the average number of particles lying at mean density in a shell of radius r and 
thickness dr inside layer /?, centered on the in-plane position of a particle in layer a. 
Quantum simulation studies of these functions by DMC methods have been reported for 
both EEB |7l8j and EHB |9l8j symmetric systems at densities ni — n2 — n and masses 
mi = ?7i2 = m, for a variety of values of the inter-layer distance d and of the coupling 
strength [here defined as — {aB^/rm)^^ where a-Q is the Bohr radius]. Theoretical 
calculations of gapir) based on a quantum STLS scheme and comparisons with DMC 
data have been reported by Moudgil et al. [8] (see references given therein to earlier 
STLS calculations). 

In the present work we develop and evaluate an alternative approach to pair corre- 
lations in bilayer systems, in the absence of inter-layer tunnel and exchange, by means 
of an extension of the analytic theory given by Asgari et al. [10 for the ground state 
of 2D quantum plasmas. This theory has provided a fully quantitative account of DMC 
data for the pair distribution function in monolayers of both electrons llj and charged 
bosons [12] , and has reproduced the location of the transition of the 2D electron gas from 
the paramagnetic to the spin-polarized state as reported in DMC work of Attaccalite et 
al. [13] . In brief, the general strategy is to set up Schrodinger-like differential equations 
for the "pair amplitudes" y/ gap{r) starting from a Fermi-hypernetted-chain (FHNC) 
approximation as first developed by Chakraborty and coworkers |14j , which is then sup- 
plemented with in-plane three-body correlation terms [15] and tailored to embody the 
Hartree-Fock (HE) limit as well as a set of sum rules for a two-component fluid of charge 
carriers. The necessary input to account for inter-layer three-body terms is regretfully 
not available for a full test of the theory at the present time. The theory is free from 
the instability that appears in the quantum STLS approach of Moudgil et al. [8J in the 
strong-coupling regime. 

The outline of the paper is briefly as follows. In Sect. [5] we present a formally exact zero- 
energy scattering equation for the pair amplitudes and introduce the approximations that 
we propose for the evaluation of the main contributions to the scattering potentials. In 
Sect. Owe pause to show that these approximations satisfy the plasmon sum rule and the 
charge neutrality condition. Our numerical results for the pair distribution functions of 
EEB and EHB systems are presented in Sect. [5] in comparison with the DMC simulation 
data ^7.9.8. . Finally, Sect. [5] summarizes our main conclusions. 
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2. THEORY 



We consider a double-quantum-well system with electrons in the first well and electrons 
or holes in the second well, with masses mi and m2 and areal densities ni and n2. We 
assume the confining potential of each layer to be so high and narrow that we can neglect 
any interlayer tunneling and treat the motion of the carriers in the two wells as strictly 
two-dimensional. The distance d between the layers is small so as to permit inter-layer 
interactions. The Hamiltonian of the system is 

Here Vaa — e^/r, Vaa = ±e^/\/f^ + <P with positive sign for EEB and negative sign for 
EHB, and (z, a) indicates the i-th particle in layer a. Our aim is to calculate the pair 
distribution functions gapif) of the bilayer, which are defined as 

5a/3W = ^-^ J2 {S{r-^(^,»)+r^,JJ•))) (2) 

O I Cry ' I'n 

^ (i.c)^U,,3) 

where the average is taken on the ground state and S is the area of the sample. 

We start the calculation, as in the work of Chakraborty and coworkers [Mj on binary 
boson mixtures and liquid metallic hydrogen, by introducing the two-component Jastrow- 
Slater variational theory involving the FHNC approximation. We make an Ansatz of the 
Jastrow-Feenberg form for the ground-state wave function as ^ = -F'ii^22^i2$i*i'2 7 where 



Pea = U^,j ^'^P 1^ ["c«a(k(,,a) - r(j,a)l)] | 



(3) 



and $Q, is a Slater determinant of single-particle plane-wave states. The HNC closure 
relations [IS] express Uajs in terms of gafs as 

Uap{r) = In gaf3{r) - [ga(3{r) - 1] + Capir) , (4) 
with Cajsir) obeying in Fourier transform the relation 

Caf3{k) = Sapik) - Saf3 - ^^[Sa-f{k) - 6a-f]Ci3-f{k) . (5) 

7 

Here Sap{k) = Sap + y/naUp /[<7q/3(?') — 1] exp (ik • y)(Py is the structure factor. Min- 
imization of the ground-state energy leads to an Euler-Lagrange equation for the pair 
distribution functions, 



-V2 + VfJr) 



\j9o.p{r) - (6) 



where the Laplacian term derives from the von Weizsacker-Herring "surface" kinetic 
energy [T7], is the reduced mass, and the effective scattering potential is given by 

Vfp{r) = vap{r) + W%{r) + Wlp{r) . (7) 

In Eq. ([7]) the scattering potential is decomposed into the sum of three terms, where the 
Bose-like term W^^ (r) contains the effect of correlations and by itself would determine 



3 



the pair distribution functions in a boson fluid, whereas the Fermi term W^^(r) derives 
from exchange and from kinetic energy terms other than the surface kinetic energy. 

The term W^^ has a very compUcated expression within the FHNC |18j. However, in 
dealing with a one-component electron fluid Kallio and Piilo [19] have proposed a simple 
and effective way to account for the antisymmetry of the fermion wave function. Their 
argument is immediately extended to a two-component Fermi fluid ^20, and leads to the 
requirement that the HF pair distribution functions are the solution of Eq. ^ in the 
high-density regime in both layers. This prescription yields in Fourier transforms 



WLik) = — / ^"/^K ^^-dh - lim W^^ik) (8) 

and W^g(fc) = 0. The rationale behind Eq. ([8]), which involves the intra-layer pair dis- 
tribution function and structure factor in the HF approximation, is as follows. The first 
term on the RHS ensures that the HF limit is correctly embodied into the theory, and the 
second term ensures that the Bose-like scattering potential is suppressed for parallel-spin 
carriers at weak coupling. 

Turning to the Bose-like scattering potential in Eq. (O, its HNC expression has been 
derived by Chakraborty [14] and its Fourier transform is given in terms of the structure 
factors of the system as 



2^n„ng 1^ ma 



/A^(fc) 



Saaik) H Saaik) 

rrin, 



(9) 

/A'(fc) 



where £a(fc) = h^k"^ / {2m a) and A{k) = 5'ii(A;)5'22(fc)-S'i2(fc)- Following Asgari et al. [10] 
we have supplemented the HNC intra-layer potential in the first line of Eq. ^ by the 
inclusion of higher-order (three-body) Jastrow-Feenberg correlations, using the theory 
developed by Apaja et al. [TS] with a weighting factor obtained by fitting DMC data on 
the ground state energy of the 2D charged-boson fluid [H] . Data for a similar inclusion 
of three-body terms in the inter-layer scattering potential are not available at present. 

It is evident that the insertion of Eqs. l!?])-® into Eq. ^ allows a self-consistent 
numerical calculation of the pair distribution functions and of the effective interactions. 
Before proceeding to the numerical solution of this problem, we examine how the approx- 
imate theory presented above fares in regard to some exact properties of pair correlations 
in bilayers. 



3. LIMITING BEHAVIORS AND PLASMON SUN RULE 

We discuss in this Section how the present theory satisfies some exact properties that 
are characteristic of charged-particle fluids. The crucial point is that the HNC expression 
for the Bose-like scattering potential exactly cancels the bare Coulomb potential at large 
distance, all other terms in the scattering potential (including the three-body correlation 
term) being instead short-ranged. This property is proved through a careful analysis of 
the effective Schrodinger equation in Eq. ([6]) by following the line of argument presented 
in detail in the work of Davoudi et al. pOj . 
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Considering for simplicity the case of a symmetric bilayer (with ni ~ n2 = n and 
nil = 1T12 = w), inversion of Eqs. ([8]) and ([9]) then yields the limiting behavior of the 
structure factors at long wavelengths as 



Siiik - 0) = ^ 0) ^ , (10) 

where fcp = ^/I-ku is the Fermi wave number, — T?]t^ j {2m\ Wpi = ^ Aime^k / m, and 
the Q's are constants to be self-consistently evaluated for each type of bilayer. The form 
of Eq. immediately ensures that the charge neutrality condition 

d^v [g^fiir) - 1] = ~5^p (11) 

is satisfied through the linear vanishing of the structure factors at long wavelengths. 

Equation (jlOp also ensures that the symmetric and antisymmetric plasmon sum rules 
are satisfied. These read 

'^*"-»>-K|(i5 '''' 

with S± = Sii± Si2 and 



c.f(fc)=.-(fc)^V^^^^i^^fc' (13) 



c.<:^fc)=<(fc) = Jl^. (14) 
V m 

In the symmetric (antisymmetric) long-wavelength mode of the EHB (the EEB) the 
carriers in the two layers move in phase (out of phase) and screen each other's potential, 
leading to acoustic dispersion relations with coefficients determined by the asymptotic 
behavior of the screened potential [21]. 



4. NUMERICAL RESULTS 



We turn to a presentation of our numerical results, which are obtained by solv- 
ing Eq. ^ through the following self-consistency cycle. We start with the trial choice 
Uapir) = g^pi''') ^iid W^p{r) = 0, and evaluate the corresponding scattering potentials 
and structure factors. A Fourier transform of the latter leads to a new value of the input 
for gap{r) and W]^p{r), and this procedure is continued till self-consistency between input 
and output is achieved. 

Figures [T] and [2] illustrate the dependence of the effective scattering potentials on the 
Coulomb coupling strength Ts and on the inverse inter-layer spacing 7 — {rsa-B,)/d for 
EHB and EEB systems, respectively. The relative phases of the oscillations shown by 
these functions in the two layers immediately indicate that inter-layer short-range order 
is growing with either of these two parameters, this state of order being in-phase for EHB 
systems and out-of-phase for EEB systems. This result, which is in full accord with the 
available DMC data [7tM5] . is illustrated in Figured 

The degree of agreement between theory and DMC data is illustrated in Figure [S] 
The inclusion of three-body correlations in the calculation of the effective intra-layer 
potentials yields almost fully quantitative agreement with the data, as was the case for 
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Fig. 1. Intra-layer (top) and inter-layer (bottom) scattering potentials V^^(?') (in Rydberg units) as 
functions of r/lrsa^) in symmetric electron-hole bilayers, for = 5 and varying inverse layer distance 
7 = {rsaB)/d (left) and for varying at 7 = 2 (right). 

the monolayer system studied by Asgari et al. [10] . In the absence of such three-body 
terms, the FHNC theory shows some quantitative discrepancies from the data in the 
inter-layer correlations, and such discrepancies are becoming severe in the EHB case as 
the two layers are brought closer together. The limit in which inter-layer excitons are 
formed in the electron-hole bilayer is evidently incorrectly described by the simple FHNC 
theory. 

5. CONCLUSIONS 

In summary, we have presented a self-consistent analytical theory of intra-layer and 
inter-layer two-body correlations in electron-hole and electron-electron bilayers, and com- 
pared its results with data from Quantum Diffusion Monte Carlo runs. The FHNC the- 
oretical approach yields quantitatively valuable results for intra-layer short-range order 
when it is supplemented with three-body correlation terms in the construction of the 
scattering potentials. There appears to be a growing need for similar three-body effects 
in the calculation of inter-layer correlations with growing inter-layer coupling. Incorpo- 
ration of the limit of inter-layer exciton formation should be sought for electron-hole 
bilayers. 
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Fig. 2. Intra-layer (top) and inter-layer (bottom) scattering potentials V^^{r) (in Rydbcrg units) as 
functions of r/{rsaB) in symmetric electron-electron bilayers, for rg = 10 and varying inverse layer 
distance 7 = (rsaB)/d (left) and for varying r-^ at 7 = 1 (right). 
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